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Abstract 


Spinel LiMn2O, offers a potentially attractive alternative to the commercialized LiCoOz 
because of its low cost, non-toxicity, good rate performance and superior safety, but 
experiences severe capacity fading as a result of Jahn-Teller cubic phase transition to unstable 
tetragonal phase especially when operated at elevated temperature (>55°C). Among the various 
techniques carried out to overcome these setbacks is doping with cations to raise the oxidation 
state of Mn thereby decreasing the amount of unstable Mn**. This study focuses on the 
investigation of the structural and electronic properties of Cu** doped LiMn2O, using ab-initio 
computational method. The structural properties of both un-doped LiMn2O4 and Cuz+ doped 
LiMn2O, were collectively obtained using both generalized gradient approximation (GGA) 
employing pseudo-potentials, plane wave basis sets and XcryDen. From the electronic 
properties of undoped LiMn2O,, it was revealed that in LiMn2O,, only Mn** ions contribute to 
conduction, even though Mn** ions and Mn** ions exist together. Increasing the number of 
Mn** brings it a unique challenge of instability. Due to this instability, avoiding the Mn** state 
of manganese prevents Jahn-teller distortion, resulting in a reduced phase transition. Based on 
the results from the first-principles calculation of Cu?* doped LiMn2O,, the Cu valence in doped 
spinel LiCuo.2sMno,7sO4 is Cu** at a fully lithiated state. When half of the Li ions are extracted, 
the Cu** ions oxidized to Cu** ions. Since Cu** is a di-valent element, the dopant act as 
electrochemically active centres for valence change, leaving two trivalent ions to only take on 
+4 oxidation states stably, thus increasing Mn-O bonding strength which ultimately suppressed 
Jahn-Teller distortions. However, previous studies show that Cu** substitution shrinks the 
lattice because the ionic radius of Cu?* is smaller than that of the host Mn** ions, thereby the 
reversible capacity of Cu-doped spinel may be reduced as Cu ions can only provide one electron 
per ion. 
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I. INTRODUCTION 


lectricity has been a necessity for humanity ever since it 

was founded, but today's global energy crisis and 
environmental pollution have made it necessary for people to 
switch from traditional energy sources to renewable ones like 
solar, wind, geothermal, and tidal energy. These renewable 
energy sources are often gathered as electrical energy before 
being distributed to consumers. However, a common issue 
with many renewable energy sources is how little electrical 
energy can be stored. Although it can be stored directly thanks 
to capacitors, the available supply is minimal compared to the 
demand for most applications [1]. Typically, electrical energy 
must be converted into another form of energy before being 
stored, therefore the choice of an energy storage technology is 
crucial. 

One of the most popular technologies that can convert 
electrical energy back into chemical energy for storage and 
continuously provide clean power for use is the use of 
rechargeable (secondary) lithium batteries. The main benefit 
of employing lithium-ion battery technology for energy 
storage is that lithium is the lightest and most electropositive 
metallic element, which makes it possible to have very high 
energy densities while being environmentally friendly [2]. 
The development of wireless and portable electronic devices 
was made possible by lithium-ion battery technology, which 
also rekindled the search for batteries to power electric and 
hybrid automobiles. 

Lithium cobalt oxide (LiCoO2) was used as the cathode in 
the initial generation of lithium-ion batteries; however, this 
substance is both expensive and hazardous. Lithium 
manganite systems appear to offer a workable solution to the 
issues as a result of research to identify an alternate material. 
A prospective replacement for currently available lithium 
rechargeable battery cathode materials based on nickel and 
cobalt oxides is the lithium manganese oxide spinel 
(LiMn20a). 

LiMn2Og, a lithium manganate with a spinel structure, has 
long been known to intercalate lithium. It has been extensively 
investigated as a positive electrode due to its exceptional 
electrochemical activity, low cost, relative cost abundance, 
environmental friendliness, stability, and high columbic 


efficiency. It also has a high potential of 3.0-4.5 V, excellent 
rate capability from its 3-D framework structure, and high 
columbic efficiency. 

However, its severe capacity fading caused by the structural 
phase transition from a stable cubic phase to an unstable 
tetragonal phase as a result of the Mn** ions’ Jahn-Teller 
distortion, disproportion reactions, and surface reactions with 
acidic contaminants in the electrolyte, especially at elevated 
temperature, severely restricts its full industrial application in 
the highly competitive energy storage market [3]. 

Doping with cations to replace a portion of Mn** ions at 16d 
sites to raise the oxidation state of Mn and thereby reduce the 
amount of unstable Mn*? is one method used to overcome 
these obstacles. Doping is an efficient way to suppress the 
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Jahn-Teller distortion and to improve the cycling performance 
of LiMn20s [4]. 

The typical method of material exploration and research 
follows the  synthesis-structure-property thread, and 
interpretations are frequently based on_ straightforward 
experimental observation and analysis. Typically, this 
procedure takes too long and is ineffective. A more effective 
and economical method has been created for forecasting and 
analysing material structures as a result of advances in 
computing technologies and _ sophisticated theoretical 
techniques using the Density Functional Theory (DFT) [5]. 

All calculations of systems at the atomic level, whether they 
are calculations to establish the atomic structure of a material 
or calculations to determine it’s electronic, optical, or other 
physical properties, begin with the fundamental laws of 
quantum mechanics. The actual unresolved problem in 
quantum mechanics is how to deal with indistinguishable, 
interacting particles, notably electrons. Density functional 
theory (DFT) is a ground-state theory that emphasizes charge 
density as the relevant physical parameter and may be 
successfully used to compute various ground state properties 
of various systems. Theoretically, the DFT can calculate any 
attribute that can be obtained by resolving Schrédinger's 
equation. It is possible to quickly screen the impact of various 
doping elements on the electrical structure of LiMn2O, using 
the first principle DFT computational modelling [6]. 

Several investigations have been conducted to improve and 

stabilize the structural characteristics of LiMn2O, by lowering 
the Mn** concentration and stabilising the cubic structure from 
the Jahn-Teller distortion effect [7]. 
By doping LiMn.O, with the Cu?* cation, this work attempts 
to identify and improve the material's structural and electrical 
properties by increasing its conductivity period. The Jahn- 
Teller phase transition is predicted to be suppressed by the 
replacement of Mn ions with Cu?* because it will increase the 
average oxidation state of Mn and reduce the number of 
unstable Mn?*. Due to its better electronic conductivity, 
capacity to decrease the activation barrier for Li diffusion in 
spinel, capacity to participate in the charge transport process, 
and ability to be more stable during electrochemical cycling, 
Cu** was chosen for this study's dopant. 


Il. COMPUTATIONAL METHOD 


The Perdew-Burke-Ernzerhof generalized gradient 
approximation (PBE-GGA) pseudo potential has been utilized 
to approximate the exchange-correlation potential for both 
doped and undoped LiMnzOzg in all computations performed 
in this study using the Quantum-ESPRESSO code. For 
integrals, the Maxfessel-Paxton smearing approach was used. 

Brillouin zone sampling (k-points) was carried out with an 
increment of 2 X 2 X 2 in the range from 2 x 2 xX 2 to10 x 
10 x 10, by fixing the other parameters constant. The values 
of minimum energies obtained to the corresponding k-point 
grids were plotted and a (6 x 6 X 6) mesh of Monkhorst-pack 
special k-points was selected for this project since the 
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structure was well converged at this k-points mesh. Also, the 
convergence test for total minimum energy as a function of 
cut-off energy (€cytwyc) is performed with an increment of 
10 Ry in the range of 20 to 60 Ry. The optimized value for 
plane wave cut-off of 40 Ry was used for all calculations. 

To ensure accuracy in this study k-point of 6 x 6 x 6 and 
cut-off energy (€cytwrc) Of 40 Ry were used to create an 
input file for the subsequent calculations. 

Doping was done by the substitution of the Mn atom by 
Cu** dopants at the (0.12500: 0.12500: 0.12500) positions. 
The four manganese atoms were replaced by one copper ion 


Total Energy 


Ecutwfc 
(a) 


Total Energ: 


making a 25% occupancy by the dopant. 


III. RESULTS AND DISCUSSION 


A. Total energy convergence test results and discussion 

The self-consistency energy convergence 
threshold (conv;y,), the Brillouin zone sampling k-points and 
the plane wave kinetic energy cut-off (Ecuewrc) for the 
material of interest are all necessary for accurate DFT 
calculations. The aforementioned parameters were optimized 
with energy reduction to acquire the right results in this work. 


Koints 


(b) 


Fig. 1 (a) convergence of the total energy of LiMn2O, to the kinetic energy cut-off (E,,,,) (b) convergence of the total energy 
of LiMn2O, to the k-point grid 


Plane-wave self-consistent field (PW,,,) calculations are 
based on convergence tests to plane-wave energy cut-off and 
k-point mesh. Fig. 1(a) illustrates the convergence of the total 
energy of LiMnzO, to the kinetic energy cut-off (E,,,¢). With 
an increment of 10Ry between20 and70Ry, the 
convergence test for total minimum energy as a function of 
cut-off energy is conducted. The remaining input file 
characteristics are fixed while the energy cut-off is varied. 

The convergence of the LiMn2Ou's total energy to the k- 
point grids is depicted in Fig. 1(b). By holding the other 
parameters constant, the convergence test for the total 
minimum energy vs. k-point sampling with an increment of 
2xX2x2x2 is carried out. The total energy F,,; (eV) 
changes significantly, as shown in Fig. 1(b), but there is some 
stability between the grid sizes (6 x 6 X 6 and 10 x 10 x 10), 
which results in less significant changes in total energy. For 
this reason, the 6 X 6 X 6 grid is chosen as the optimal k-point 
value. 
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B. Un-doped Structural and Electronic Properties 


Fig. 2 shows that the stable phase of undoped LiMn204 
crystallizes in the cubic Fd3m space group at room 
temperature. In this cubic structure, manganese ions occupy 
the octahedral positions (16d), oxygen ions occupy the 
tetrahedral positions (32e), and lithium ions occupy the 
tetrahedral (8a) positions. Each of them has voids separating 
it from four nearby ions (16c). The remaining vacant spaces 
are connected and provide a 3D diffusion channel for the quick 
diffusion of lithium ions. The three-dimensional MOz host 
serves as the structure's foundation, and the transition metal 
layer's vacancies guarantee the three-dimensional Li diffusion 
paths. 

Another finding in the structure mentioned above is that all 
of the atoms in the spinel LiMn2Os structure are linked by 
strong chemical bonds that run across the lattice, making it 
simple for Li to intercalate or de-intercalate into the structure. 
Each 8a tetrahedron shares all four of its faces with a l6c 
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octahedral unit, which is a significant characteristic of this 
cubic spinel structure. 


Fig. 2 Structure of Un-doped LiMn204 using XcryDen. 


C. Density of state (DOS) of Un-doped LiMn204. 

The electronic structure of undoped LiMn2O, is examined 
in this section using the density of state (DOS). Fig. 3 and 4 
display, respectively, the DOS of LiMnzO,, Li ions, Mn ions, 
and O ions. While the lithium atoms in the spinel structure are 
significantly ionized, the electronic structure of LiMn2Oy, is 
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controlled by a strong hybridization between the Mn-d and the 
O-p states. 

With an equal amount of anti-ferromagnetic aligned Mn** 
and Mn‘ ions, the average valence of Mn ions is +3.5. The d 
orbitals split into two E, orbitals as a result of the Mn ions in 
the O octahedral sites. Because there is a band gap 
betweent 2, and e, and three valence electrons filling three 
tz, orbitals while leaving one unfilled in Mn", the ion acts as 
an insulator. 

According to Fig. 4, the 3d states of Mn atoms contributed 
to the first peak in the valence band, which is dominated by 
these states along with some contributions from the 3p states 
of O atoms and the 2s states of Li. There are no maximum 
peaks in the conduction band; instead, the shallow peaks just 
before the Fermi energy are dominant. These peaks come from 
the 3p orbitals of oxygen and are dominated by three broad 
peaks. 

Thus, from this research, even though Mn** ions and Mn** 
ions coexist in LiMn2O4, only Mn3* ions contribute to 
conduction. But as the Mn** content rises, instability poses a 
special problem. This instability causes Jahn-teller distortion 
to be prevented by the manganese in the Mn** state, which 
reduces phase transition. 
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Fig. 3 The density of state of Un-doped LiMn20,4 
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Fig. 4 Partial Density of State of Un-doped LiMn2O, indicating the number of states per atom. 


D. Band Structure of Un-doped LiMn204 

The calculated electronic band structures of Un-doped 
LiMnzOg are shown in Fig. 5, using pbe-GGA, the bands were 
selected along 6 symmetry points which are (IT > X — K > 
L — X —T) and the band structure energy range is plotted 
from 0.0 eV to 15.0 eV. The red line represents Fermi Level. 
LiMn204 is considered a semiconductor based on the above 
calculation of its band structure. Also, as seen from Fig. 5, at 
I point, the presence of the valence band and conduction band 
around different symmetry points. This indicates that the 
materials have indirect band-gap nature. 

The resulting electronic band structure resulting from the 
calculation exhibits a small band gap of 4.3 eV, in agreement 
with previous studies using similar parameters [8]. 

15 


40F 


E-Et(eV) 


Fig. 5 Calculated band structures of the undoped LiMn204 
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E. Structural and electronic properties of Cu?” Doped 
LiMn0,4 

The schematic primitive diagram of the Cu’* doped 
LiMn20, lattice atom group is shown in Fig. 6. The face- 
centre cubic lattice makes up each of the two lithium, three 
manganese, one copper, and eight oxygen atoms that make up 
the structure of Cu2* doped LiMn2Ou. 

The cubic spinel LiMn20s has a cation distribution of (Li)8a 
(Mnz)16d O4 as seen in Fig. 6, the subscripts 8a and 16d 
denote Fd3m symmetric cation sites. Tetrahedral 8a sites are 
occupied by Li ions, while the octahedral 16d sites are 
occupied by Mn ions. The O ions at the 32e sites form a cubic 
closest packed (ccp) structure. 


Fig. 6 Structure of Cu2+ doped LiMnzO, using XcryDen. 
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Mn ions have an average valence state of +3.5, whereas Li 
and O ions have valence levels of +1 and -2, respectively. At 
the 16d octahedral location, some of the Mn ions are swapped 
out for Cu* ions, creating a solid solution with the chemical 
formula LiCuo.2sMnq.75Os. 


F. Density of state (DOS) of Cu?* doped LiMn2O4 

The Density of State of Cu ** doped LiMn2O, is shown in 
Fig. 7 the highest peak in the valence band denotes 
contributions from the Mn atom's 3d states. 

In the conduction band, the Mn-4d orbitals form the small, 
high peak immediately after the Fermi energy, and the visible 
shallow peaks are states derived from a combination of Mn- 
3d, Cu-1s, and O-2p orbitals. The sharp peaks just before the 
Fermi energy originates from the Cu-3d orbitals. 

Fig. 8 and 9, which exhibit the partial density of state of the 
dope LiCuo.2sMn;.75O4 and the partial density of state of the 
doped LiCuo.soMn;.s004 respectively, demonstrate the number 
of states per atom. 

Li* ions, Mn** ions, and Mn** ions are the three main forms 
of cations found in LiMn2O.. In Fig. 8, a 25 per cent Cu doping 
was performed. It can be seen that the first and highest peak 
results from a Mn-|s state, followed by a Mn-3p state and a 
Cu-3d state just before the Fermi energy, and the conduction 
band. From the observation that Mn-3p results in the highest 
peak at the conduction band, it can be inferred that Cu ions 
occupy 4b octahedral sites and Li ions. 
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The highest peak in Fig. 9, which has more Cu doping, is 
caused by the Mn-3p state; in contrast to Fig. 7, which clearly 
shows that there are more Mn** ions than Mn** ions, Li ions 
are more likely to be surrounded by Mn** ions. Due to the 
electrostatic effect, combining Li* and Mn** ions can limit the 
localization of positive charges and further lower the system's 
overall energy. 

As aresult, Mn** ions have less dense electron clouds than 
Mn** ions, which results in weaker Mn-O interactions and 
longer Mn-O bond lengths. These shorter Li-O bond lengths 
then result in smaller Li 16c site octahedral volumes. Due to 
the stronger Coulombic repulsion between Li*-Mn** and Li*- 
Mn**, the Li*-Mn* distances are longer than Li*-Mn** 
distances [9]. When the electrostatic balance is reached, the 
Li-Mn**/Mn** distances spread with wide distribution, 
reducing the elongation or contraction of the Mn-O bond 
pointing towards vertex O as a result of the asymmetric shape 
of d electron 

These first-principles calculations' findings indicate that the 
Cu valence in the doped spinel LiCuo.2sMno,.7sO4 is Cu2* in the 
fully lithiated state. The Cu** ions oxidized to Cu** ions upon 
extraction of half of the Li ions. However, because Cu ions 
can only supply one electron per ion, the reversible capacity 
of Cu-doped spinel may be diminished. But it is possible to 
lower the Li diffusion activation barriers when doped with a 
small amount of Cu, as in Fig. 8, to achieve a higher rate 
capability. 
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Fig. 7 Density of State of Cu* dope LiMn.O4 
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Fig. 8 Partial Density of State of dope LiCuo.2sMnj.75O4 indicating the number of states per atom 
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Fig. 9 Partial Density of State of dope LiCuo.soMnj.s0O4 indicating the number of states per atom 


G. Band Structure of doped LiMn204 

Fig. 10 shows the calculated electronic band structures of 
doped LiCuo2sMn;7sO4 the bands were selected along 6 
symmetry points which are (T — X — K—> L— —T) and 
the band structure energy range is plotted from 0.0 eV 
to 15.0 eV. The red line represents Fermi Level. Based on the 
above calculation of band structure, LiMn 2O, is considered a 
semiconductor. Also, it shows the presence of the valence 
band and conduction band at around different symmetry 


VOLUME 03, ISSUE 01, 2023 


points which indicates that the materials have indirect band 
gap nature. 

The estimated results are consistent with some experimental 
data since the resulting electronic band structure has a wide 
valence bands width (VB) of 6.90 eV and a small conduction 
bandwidth (CB) of 2.24 eV. The electronic band structure 
calculated here is fairly similar to the reported values of 
1.923 eV [10]. 
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Cu** is a divalent element, so the dopant functions as 
electrochemically active centres for valence change, leaving 
only two trivalent ions to take on stable +4 oxidation states. 
This increased the strength of the Mn-O bond, which in turn 
reduced Jahn-Teller distortions. Cu** substitution, however, 
causes the lattice to contract because its ionic radius is less 
than that of the host Mn3* ions [4]. 


E-Ef (eV) 


Fig. 10 Calculated band structures of the doped LiMn204 


IV. CONCLUSION 


Numerous significant electronic parameters, including band 
gaps and density of states, have been measured for both 
undoped and Cu** doped LiMn20, and results compared with 
theoretical and experimental findings. Though Mn** ions and 
Mn** ions coexist in LiMn20Oy, it was discovered from the 
electrical characteristics of undoped LiMn2Oy, that only Mn?* 
ions contribute to conduction in this material. However as the 
Mn** content rises, instability poses a special problem. Due to 
its instability, avoiding the manganese's Mn** state is 
necessary to avoid Jahn-Teller distortion and a shortened 
phase transition. The Cu valence in doped spinel 
LiCupo,2sMnp.75O4 is Cu** ina fully lithiated state, according to 
the findings of a first-principles calculation of Cu’* doped 
LiMn2O,. After removing half of the Li ions, The Cu** ions 
oxidized to produce Cu** ions. Cu** is a divalent element, so 
the dopant functions as electrochemically active centres for 
valence change, leaving only two trivalent ions to take on 
stable +4 oxidation states. This increased the strength of the 
Mn-O bond, which in turn reduced Jahn Teller distortions. 
However, Cu?" substitution causes the lattice to contract 
because the ionic radius of Cu?* is smaller than that of the host 
Mn** ions. As Cu ions can only provide one electron per ion, 
this may reduce the reversible capacity of Cu-doped spinel. 
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